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1. Statement 

Novelty (N) 

Inventive step (IS) 
Industrial applicability (lA) 



Claims 
Claims 

Claims 
Claims 

Claims 
Claims 



2-9. 11-17 
1.10,18 



1-18 



1 - 18 - 
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NO 

YES 
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2. Citations and e>q)lanations (Rule 70.7) 

The problem solved is a more efficient way of calculation in 
converting measurement signals emitted from electromagnetic 
fields . 

Reference is made to the following documents: 
D1:US5269325 
D2:EP10302160 
D3:US6370414 
D4:US63 74131 
D5:EP0982597 

D6:EP0632353 . ^. ^, . ^ 

D7- -Multipole expansions of electromagnetxc fxelds usxng 
Debye potentials', C. G. Gray, American Journal of 
Physics, Vol. 46, pp. 169-179, 1978. 



Dl through D6 each disclose the method and the device of 
determining current distribution of an object by measuring the 
magnetic fields in the vicinity of the object, subsequently 
calculating the signals as from a set of virtual sensors. 

D7 discloses a method of multipole expansions of 
electromagnetic fields. 

The invention according to claims 1, 10 and 18 is stated by 
the technical features: (i) converting signals into orthogonal 
signals from a set of virtual sensors, and (ii) determining 
the current distribution based on the virtual sensors. All 
technical features are known from each of Dl through D6. The 
applicant states that D5 and D6 are irrelevant. However, D5 
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discloses a method of measuring weak signals and applying 
virtual channels [0044] . D6 describes the measurement of a non 
biological object. Present claims 1, 10 and 18 are not 
directed to any particular kind of object. D6 further 
describes virtual sensors. Accordingly, all cxted documents 
are relevcuit. 

The applicant also states that the source model needs no 
further regular isat ion and that no statistical analysis is 
required. However, none of the technical features m the 
independent claims implement those effects. 

AS to the argument that the equations in D7 are not presented 
in a straightforward maimer, this is neither the case in the 
independent claims. On the contrary, claim 1 only contains the 
desideratum of «a predetermined function basis to be 
efficiently calculated" . 

The invention claimed in claims 1, 10 and 18 does consequently 
not fulfil the requirement of novelty. 

The invention according to claims 2 and 11 differs from what 
is known from for instance Dl, regarded as the docxament that 
is closest in describing the invention, by stating the use of 
raultipole expansion of the field. The problem solved is an 
alternative way of calculation that is faster and makes the 
mathematical operations simpler. 

Paced with the problem of an alternative way of analyzing the 
signals, the person skilled in the art knows a solution from 
D7. D7 addresses multipole expansions of electromagnetic 
fields. Since Dl and D7 are both applied to electromagnetic 
fields! it is considered obvious that a person skilled in the 
art would combine the information and arrive at the invention 
claimed in claims 2 and 11. 

The features stated in the remaining claims are technical 
features considered known from the cited documents or are 
general common knowledge for a person skilled in the art, for 
instance claims 3-6 and 12-14. Claims 7-9 and 15-17 state 
features that differ from what is disclosed in D7. However, 
the difference does not exhibit a technical character. The 





Intemational application No. 1 


INTERNATIONAL PRiaLIMINARY REPORT ON PATENTABILITY 


PCT/FI2004/000532 | 


Snpplemental Box 


In case the space in any of the preceduig boxes is not sufficient. 




Continuation of: Box V 





difference is merely a mathematical theory, and as such it 
does not involve an inventive step (PCT Rule 39 (x) ) . 

Accordingly, the inventions claimed in claims 2-9 and 11-17 do 
not fulfil the requirement of inventive step. 



The invention is industrially applicable. 



g^saSo^agj^r] ?ct / fi 2oow o o o 5 

1 2 -12- 200s 

1 

METHOD FOR USING A MOLTI-CHANNEL MEASUREMENT SXGNAI. IN 
SOURCE MODSLI1IN6 

FIELD OF THE INVENTION 

The present invention relates to a new and 
advanced method for using a multi-channel measurement 
signal in source modelling. In particular^ the present 
invention relates to a new way of converting the meas- 
urement signals measured using a multi-channel measur- 
ing device that measures a curl free and divergence 
free vector field into a form optimal from the stand- 
point of the source modelling. 

BACKGROX7ND OF THE INVENTION 

The electrical operation of an object to be 
examined can be examined e.g. by measuring the mag- 
netic field produced by the electric currents of the 
object using sensors placed outside the object. The 
modelling of a source distribution thus obtained based 
on the measured signals is^ however, very difficult 
because each magnetic field distribution can be pro- 
duced by many different source distributions. In other 
words, a source distribution cannot be solved unambi- 
guously based on the measured signals, so to solve the 
problem, different restrictive conditions must be set, 
such as some parametric model based on prior informa- 
tion for a current, or a non-parametrised norm re- 
striction. 

For the non-parametrised modelling of a con- 
tinuous current distribution, minimum norm estimates 
are usually used, in which there is an attempt to ex- 
plain a signal measured using a multi-channel measure- 
ment device with a current distribution whose norm is 
as small as possible. As the norm, usually an LI or L2 
norm is selected, of which the previous is the sum of 
the lengths of the current elements over the selected 
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volume, and the latter is the sum. of the squares of 
the lengths of the current elements over the selected 
volume. The calculation of the minimtam norm estimates 
has been described e.g. in publications "Interpreting 
magnetic fields of the brain: minimum norm estimates", 
M.S. Hamaiainen et al. Medical & Biological Engineer- 
ing & Computing, Vol. 32, pp. 35-42, 1994, as well as 
"Visualization of magnetoencephalographic data using 
minimum current estimates", Uutela K. et al, Neurol- 
mage. Vol. 10, pp. 173-180, 1999. 

Conventional minimum norm estimates involve 
inherent problems such as slowness of calculation and 
susceptibility to noise. For example, in the case of 
an L2 norm, one needs an inverse matrix of matrix G, 
whose element (i, j) contains the inner product of the 
lead fields of the ith and jth measurement sensor, so 
one must calculate these inner products for each pair 
of sensors. The lead field is so determined that the 
signal measured by a sensor is the projection of the 
current distribution for the lead field of the sensor 
in question. The noise problems are due to the fact 
that matrix G calculated for the sensors is suscepti- 
ble to noise, so in the calculation of its inverse ma- 
trix, regularisation is needed in the practical situa- 
tions . 

Regularisation methods, such as the 
truncation regularisation of the singular value decom- 
position, usually are non-intuitive, and usually also 
to be solved for each case specifically. A regularisa- 
tion of the wrong type may lead to an erroneous model- 
ling result. 

Therefore, source modelling nowadays still 
involve problems such as the hardness and slowness of 
the computation, the possible errors caused by noise, 
as well as the case-specificity due to the regularisa- 
tion. Further, as stated above, the regularisation may 
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cause considerably errors to the final computation re- 
sult. 

OBJECTIVE OP THE INVENTION 

The objective of the invention is to elimi- 
nate the drawbacks referred to above, or at least sig- 
nificantly to alleviate them. One specific objective 
of the invention is to disclose a new type of method 
which can be used to considerably lighten and acceler- 
ate the computation associated with the modelling of a 
continuous current distribution, as well as to lessen 
the problems with noise. 

As for the features characteristic of the in- 
vention, reference is made to them in the claims. 

DESCRIPTION OP THE INVENTION 

The present invention relates to a new kind 
of manner of . determining the continuous current dis- 
tribution of an object being examined using basis vec- 
tor components of signal space calculated from the 
measured signals. The components in question have been 
so selected that they describe the features, as inde- 
pendent as possible, of the current distribution being 
examined, which enhances the computation and makes it 
more accurate. 

The basic idea of the invention is that be- 
cause the computation of the inner products of the 
sensor fields is heavy and difficult using a conven- 
tional set of sensors, it is worth using a special set 
of sensors, whose lead fields are orthogonal and, if 
possible, to be analytically computed. In principle, 
this can be implemented by a suitable physical set of 
sensors. As a suitable physical set of sensors is, 
however, often quite difficult to manufacture, it is, 
in most of the cases, advantageous to use virtual sen- 
sors computationally generated from a conventional set 
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of sensors, i.e. the measurement signals are converted 
into other ones by a suitable conversion so that they 
correspond to the signals that the virtual measurement 
device would have measured. At the same time, it is 
possible, if necessary, to eliminate the signals asso- 
ciated with external interference. This conversion has 
been described e.g. in patent application FI20030392, 
which is incorporated herein by reference. After the 
conversion, the source modelling is performed in an 
optimal manner using the basis vector components of 
the signal space instead of the actual measurement 
signals. One substantial feature of the invention is 
that after the conversion, the source model need not 
be any more regularised. 

Thus, the present invention relates to a 
method for determining the current distribution of an 
object by measuring the. magnetic fields in the vicin- 
ity of the object using a multi-channel measurement 
device. Advantageously, at least one measurement sen- 
sor corresponds to each channel, and the object is ap- 
proximated by a spherically symmetric conductor. The 
object can be e.g. a human being's head. 

According to the invention, a multi-channel 
measurement signal corresponding to each measurement 
sensor is converted into the signals of a predeter- 
mined set of virtual sensors, and the current distri- 
bution of the object being examined is determined by 
depth r from the signals of the set of virtual sensors 
in a beforehand selected orthonormal function basis. 
In that case, the estimation of a current distribution 
is fast and robust. Further, to achieve the set of 
signals corresponding to the set of virtual sensors, a 
multi-pole expansion is calculated from a multi- 
channel measurement signal. A multi-pole expansion can 
be calculated in two ways: by taking into account the 
magnetic fields emitted by sources outside the object 
being measured, or by ignoring them. 
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Advantageously, as the orthonormal function 
basis, a basis with the following fo3nn is selected: 



/.O m— / 



wherein fi(r) is a radial function to be 
freely . selected and X,^(.0>(p) is so-called vector 
spherical harmonic. In that case, it is possible to 
place the function basis into a current distribution 
equation, and the coefficients of the current distri- 
bution are analytically solved based on the equation: 



wherein is a constant associated with or- 
der 1 and R is the radius of the sphere to be exam- 
ined. Advantageously, function fx(r) is used to adjust 
the depth weighing of a current distribution model. 

Furthermore, the invention relates to a meas- 
urement device for determining the current distribu- 
tion by measuring magnetic fields in the vicinity of 
the object. The measurement device includes a set of 
measurement channels that measure a curl free and 
divergence free vector field, whereby at least one 
measurement sensor corresponds to each channel, and 
processing means for processing the measurement sig- 
nal. Advantageously, the object is approximated with a 
spherically symmetric conductor. 

According to the invention, the processing 
means include a conversion module for converting a 
multi-channel measurement signal corresponding to each 
measurement sensor into the signal of a predetermined 
set of virtual sensors; and calculation means for de- 
termining the current distribution of an object being 
examined or for calculating by depth r from the sig- 
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nals of the set of virtual sensors in a beforehand se- 
lected orthonormal function basis. In one embodiment, 
the calculation means are arranged to calculate a 
multi-pole expansion from a multi-channel measurement 
signal. 

The invention enables one to considerably 
lighten and accelerate the calculation associated with 
the modelling of a continuous current distribution. 
The invention further enables one to reduce the prob-r 
lems associated with noise. Further, the invention 
simplifies the regularisation of a source model, or 
eliminates the need for it, thus also significantly 
reducing the possibility of error. 

LIST OF FIGURES 

In the following, the invention will be de- 
scribed in detail by means of examples with reference 
to the accompanying drawing, in which 

Fig. 1 represents one measurement device in 
accordance with the present invention; and 

Fig. 2 is a flow chart illustrating one 
method in accordance with the present invention. 

DETAILED DESCRIPTION OF THE INVENTION 

Fig. 1 represents one measurement device in 
accordance with the present invention comprising a set 
of measurement channels iS iS . . . 1^ that measure a 
curl free and divergence free vector field, in which 

1. 2 5 

case at least one measurement sensor 2 , 2 ,...2 cor- 
responds to each channel; and processing means 3 for 
processing the measurement signal. Advantageously, the 
processing means have been implemented by means of a 
computer. Further, the processing means include a con- 
version module 4 for converting a multi-channel meas- 
urement signal corresponding to each measurement sen- 
sor into the signals of a predetermined set of virtual 
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sensors, and calculation means 5 for determining the 
current distribution of an object being examined. 

Further Fig. 1 shows the object being meas- 
ured K, into the vicinity of which the measurement 
sensors 2^ 2^ , . . .2^ have been placed. Inside the ob- 
ject being measured there is presented a current loop 
which describes the source of the magnetic field being 
measured. The object being measured can be e.g. a hu- 
man head, and there can be several current sources. 

Fig. 2 represents the main steps of one pre- 
ferred embodiment of the present invention. At first, 
a set of sensors is used to measure a multi-channel 
measurement signal, step 21. After that, the signal is 
converted to correspond to the signal measured using a 
so-called set of virtual sensors, step 22, making the 
mathematical operations simpler. Finally, it is possi- 
ble to simply calculate the current distribution in an 
object from the signal of the set of virtual sensors, 
step 23, i.e. in practice, to describe the places and 
strengths of the current loops inside a sphere or 
head. 

In the following section, the mathematical 
background and grounds of the invention are described. 
When the magnetic fields are converted into coeffi- 
cients M,„='a^+ib,„ associated with the basic solution 
f.<M)Ylm(0,<p) of the Laplace equation, wherein i is an 
imaginary unit, they can be expressed by means of the 
current distribution J(f) in spherical coordinates 
{r,0j(p) , whereby they are of the form: 

Afto =r, yX^{e,<p).J(f)dv (1) , 

V 

wherein the integration is performed over the 
volume being examined^ Yi is a constant associated 
with order 1 and Xlmi9,q>) is so-called vector spherical 
harmonic. 
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This form can be derived e.g. from publica- 
tion "Multipole expansions of electromagnetic fields 
using Debye potentials", C.G. Gray, American Journal 
of Physics, Vol. 46, pp. 169-179, 1978. The expression 
mentioned above is of the lead field form, wherein the 
lead field of the multi-pole coefficient Mm is of the 

form: 

4,(0 = ^'^/«(^^<^) (2) . 

On the other hand, the vector spherical har- 
monics form by depth r an orthonormal basis, so with 
the depth in question, the current distribution can be 
presented in the function basis in question: 

= S Z^/«// ir)Xlmi0,9>) ( 3 ) , 

wherein fi(r) is some radical function. 

When as the volume to be examined, spherical 
volume is selected, by placing the previous expression 
into the equation (1), the coefficients of the current 
distribution can be analytically solved: 

wherein is a constant associated with or- 
der 1 and R is the radius of the sphere to be exam- 
ined. The previous equation (4) indicates that the co- 
efficients of a current distribution model presented 
in an orthonormal basis can be solved based on coeffi- 
cients Mim in a completely trivial manner using ana- 
lytical expressions without any kind of regularisa- 
tion. This is computationally very fast and nvimeri- 
cally stabile. Function fi(r) is freely selectable. 
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and can be used to adjust the depth weighing of a cur- 
rent distribution model. 

Furthermore, it must be noted that in the 
case of a spherical conductor, it is possible to use 
as basis functions also some other orthogonal basis or 
a basis whose inner products can be otherwise quickly 
calculated. This is achieved e.g. by slightly breaking 
the orthogonality in a manner known per se, or by us- 
ing a basis which is not orthogonal but whose inner 
products can be calculated beforehand. 

The invention is not limited merely to em- 
bodiment examples referred to above, instead many 
variations are possible within the scope of the inven- 
tive idea defined by the claims. 



AMENDEDSHEET 



The Swedish Patsnt Office 
POT Matematton^ Appllcatio a 



^CT /Fl 200^}/ 00 0 53 
10 1 2 -12- 2005 



CIi2\lMS 



1. A method for determining the current dis- 
tribution of an object by measuring the magnetic 
fields in the vicinity of the object using a multi- 
channel measurement device that measures a curl free 
and divergence free vector fields whereby one measure- 
ment sensor corresponds to each channel, charac- 
terised in that 

converting a multi-channel measurement signal 
corresponding to each measurement sensor into the sig- 
nals of a predetermined set of virtual sensors, which 
signals are mutually orthogonal; and 

determining the current distribution of an 
object being measured from the signals of the set of 
virtual sensors in a predetermined function basis to 
be efficiently calculated. 

2. The method as defined in claim 1, 
characterised in that the object is- ap- 
proximated using a conductor, and a multi-pole 
expansion of the field is calculated from the multi- 
channel measurement signal, 

3. The method as defined in claim 2, 
characterised in that the object is ap- 
proximated using a spherically symmetric conductor. 

4. The method as defined in claim 2, 
characterised in that the multi-pole 
expansion is calculated by taking into account the 
magnetic fields emitted by sources outside the object. 

5. The method as defined in claim 2, 
characterised in that the multi-pole 
expansion is calculated by ignoring the magnetic 
fields emitted by sources outside the object. 

6. The method as defined in claim 2, 
characterised in that the external inter- 
ferences are eliminated using some other known inter- 
ference eliminating method prior to the conversion. 
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7. The method as defined in claim 2, 
characterised in that as the orthonormal 
function basis, a current distribution equation of the 
following form is selected: 

wherein cim are coefficients of the current 
distribution, fx (r) is a freely selectable radial 
function and Xlm{0,<p) is so-called vector spherical 
harmonic- 

8. The method as defined in claim 5, 
characterised in that 

the orthonormal function basis is placed into 
a current distribution equation; and 

the coefficients of the current distribution 
are analytically solved from the equation: 

wherein is a constant associated with or- 
der 1, Mi^ are multi-pole coefficients, R is the radius 
of the sphere to be examined, and fi(r) is a freely 
selectable radial function. 

9. The method as defined in claim 7, 
characterised in that function fx (r) is 
used to adjust the depth weighing of the current dis- 
tribution model. 

10. A measurement device for determining the 
current distribution of an object by measuring mag- 
netic fields in the vicinity of the object, the meas- 
urement device comprising: 

a set of measurement channels (1, l^, ^ 
. . . l*") that measure a curl free and divergence free 
vector field, in which case at least one measurement 
sensor 2, 2^, 2^, ...2^ corresponds to each channel; and 
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processing means (3) for processing the meas- 
urement signal and in which the object is approximated 
.using a conductor, characterised in tha: 

the processing means include a conversion 
module (4) for converting a multi-channel measurement 
signal corresponding to each measurement sensor into 
the signals of a predetermined set of virtual sensors, 
which signals are mutually orthogonal; and 

calculation means (5) for determining the 
current distribution of an object being examined from 
the set of virtual sensors using depth r in a prede- 
teinnined orthonormal function basis. 

11. The measurement device as defined in 
claim 10, characterised in that the calcu- 
lation means (5) are arranged to calculate a multi- 
pole expansion from the multi-channel measurement sig- 
nal. 

12. The measurement device as defined in 
claim 10, characterised in that the object 
is approximated using a spherically symmetric conduc- 
tor. 

13. The measurement device as defined in 
claim 11, characterised in that the multi- 
pole expansion is calculated by taking into account 
the magnetic fields emitted by sources outside the ob- 
ject being measured. 

14. The measurement device as defined in 
claim 11, characterised in that the multi- 
pole expansion is calculated by ignoring the magnetic 
fields emitted by sources outside the object Joeing 
measured. 

15. The measurement device as defined in 
claim 11, characterised in that as the or- 
thonormal function basis, a current distribution equa- 
tion with the following form is selected: 
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wherein ci„ are coefficients of the current 
distribution, fi Cr; is a radial function to be freely 
selected and Xlm(0,<p) is so-called vector spherical 
harmonic. 

16. The measurement device as defined in 
claim 14, characterised in that 

the orthonormal function basis is placed into 
the current distribution equation; and 

the coefficients of the current distributi 
are solved analytically from the equation: 



on 



wherein is a constant associated with or- 
der 1, Mi„, are multi-pole coefficients, R is the radius 
of the sphere to be examined and fj.(r) is a radial 
function to be freely selected. 

17. The measurement device as defined in 
claim 15, characterised in that function 
fi(r) is used to adjust the depth weighing of a cur- 
rent distribution model. 

18. The measurement device and analysis soft- 
ware as defined in claim 10, wherein the measurement 
device converts the signals into a set of virtual sen- 
sors prior to their storage, and the analysis software 
converts the stored data into a current distribution. 
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